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I review the latest heavy flavor measurements at RHIC experiments. 
Measurements from RHIC together with preliminary results from LHC offer 
us an opportunity to systematically study the sQGP medium properties. 
In the end, I will outlook a prospective future on precision heavy flavor 
measurements with detector upgrades at RHIC. 

PACS numbers: 25.75.-q 

1. Introduction 

Heavy quarks are expected to be a clean and penetrating probe to study 
the sQGP matter created in heavy ion collisions because of its intrinsic 
large mass property. By studying interactions between heavy quarks and 
medium, one can learn in detail the flavor dependence of parton energy loss 
mechanism, and also the medium's degree of thermalization by looking at 
the medium response to heavy quarks. Precision measurements on heavy 
flavor hadrons in a wide kinematic region will be unique to understand these 
details. Full reconstruction of charm hadrons has significant advantages 
over semi-leptonic decay leptons because of complete kinematics and clean 
interpretations. 

Quarkonium suppression due to color screening has been originally pro- 
posed as a smoking gun for QGP formation. But various cold and hot 
medium effects complex this story and there hasn't been a direct evidence 
of color screening so far. Bottomonium (or T) production at RHIC may of- 
fer us a unique opportunity to directly observe this signature because many 
cold and hot medium effects (cold nuclear absorption, regeneration etc.) are 
expected to negligible at RHIC. While at LHC, these become complicated, 
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particularly the regeneration process for bottomonium production can be 
significant, which introduces difficulties in the interpretation. 

2. Latest RHIC Results 

In this proceedings, I would like to focus on latest RHIC results on charm 
hadron and quarkonium (J/V', measurements. 

2.1. Charm hadron measurements 

Figure [l] shows the recent STAR measurements on the charm hadron 
production cross sections in p + p and minimum bias Au+ Au collisions [1] . 
Charm hadrons were reconstructed via hadronic decays. Although it was 
not possible to reconstruct secondary decay vertices for charm hadrons in 
these measurements, STAR managed to overcome the large combinatorial 
background with large amount of statistics. The left plot shows the mid- 
rapidity production cross sections for and D*^ scaled to cc pairs vs. pT 
and the result is compared to a pQCD FONLL calculation [2]. One can see 
over a wide px region, the measurement is consistent with the upper bound 
of this FONLL calculation. An interesting finding is that the charm hadron 
cross sections measured by CDF [3j and ALICE [4J at higher energies up to 
7 TeV are also closer to the upper limits of FONLL calculations. Similarly, 
the cross section of non-photonic electrons is also consistent with the upper 
bound of FONLL calculations at pT{e) > 1 GeV/c [HE]. The STAR D 
meson measurement covers about 70% of total pT acceptance, leading to 
a reasonable constrain to the total charm cross section. Fig. [T] right plot 
shows the total charm cross section per nucleon-nucleon collisions at mid- 
rapidity from p + p to central Au-|-Au collisions. Latest STAR Au-|-Au 
results were extracted from the spectrum measurements covering pT 
from 0.4 up to ~ 5 GeV/c. The cc cross sections were obtained assuming 
the c — > fragmentation ratio still holds in Au-|-Au collisions which need to 
be tested in future measurements. The results exhibit an approximate Nf,in 
scaling indicating charm quarks are predominantly produced from initial 
hard scatterings. 

The spectrum in Au-|-Au collisions can be compared to the reference 
p + p data to calculate the nuclear modification factor Raa- Figure [2] left 
panel shows the recent results from the STAR measurements. The Raa 
covers px up to ~ 5 GeV/c. Although the values are consistent with unity 
given current uncertainties, there seems to be a hint of modification for D'^ 
production in this pr region. The blue dashed line depicts the Raa calcu- 
lated based on the Blast- Wave fit to the Au-|-Au spectrum. If one uses the 
Blast- Wave model with light hadron freeze-out parameters and calculates 
the Raa, the result is shown as shaded area in the plot. The data points 
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Fig. 1. Left: Charm hadron [D'^ and production cross section inp+p collisions 
at y/s ~ 200 GeV from STAR. Measurements are consistent with the upper bound 
of the FONLL pQCD calculation. Right: Total charm production at mid-rapidity 
per nucleon-nucleon collisions from p-\- p to central Au-I- Au collisions from STAR. 
The measurements demonstrate approximate number of binary scaling for the total 
charm production cross section. 



are significantly different from this predication indicating the charm hadrons 
freeze-out differently from the system compared to light hadrons. When go- 
ing towards higher p-p, the data points tends to show a slight suppression, 
however, the uncertainties are large. The right panel of Fig. [2] plots the 
Raa from STAR together with the high pT measurement from central 
(0-20%) Pb-hPb coUisions by ALICE [4J. One should note the difference in 
the centrality for these two sets of data points. After considering this, there 
seems to be a consistent trend of Raa between RHIC and LHC although 
the uncertainties need to be shrink for detailed investigation. One should 
be aware that Raa less than unity doesn't necessarily mean the suppres- 
sion of charm quark production in heavy ion collisions. The sQGP medium 
in heavy ion collision may modify the distributions of charm quarks into 
different charm hadrons compared to p + p collisions because hadronization 
scheme other than fragmentation (e.g. coalescence) can be significant as 
already been observed for light flavor hadrons. A complete understanding 
of charm quark energy loss requires measurements of all ground state charm 
hadrons in heavy ion collisions. 



2.2. J/ij) measurements 

J /ip production has been reported by the PHENIX collaboration, focus- 
ing on the low pT suppression observation [7] . With a couple of subsystem 
upgrades, STAR has been able to measure J/ij^ production with improved 
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Fig. 2. Left: meson Raa in 0-80% minimum bias Au-|-Au collisions at RHIC 
from STAR. The blue dashed line depicts the Raa with the Au-f Au Z?" spectrum fit 
to the Blast- Wave model. The shaded area depicts the Raa with Au+Au spectrum 
calculated in the Blast- Wave model with light hadron freeze-out parameters. Right: 
D° meson Raa in 0-80% Au-I-Au collisions at = 200 GeV fr-om STAR and 

D meson Raa in 0-20% Pb-^Pb collisions at y^J]^ = 2.76 TeV from ALICE. 



statistics, particularly with significant capability covering high pj- up to 
~ 10 GeV/c. A combination of botfi PHENIX and STAR measurements 
covers significant large acceptance in 47r [HI El [TU], offering us a great op- 
portunity to constrain the quarkonium production mechanism as well as to 
learn the cold and hot medium properties. Most newly developed models 
can produce the cross section data points within accessible kinematic region. 
Precision cross section measurements provide constrains on model calcula- 
tions. In the meantime, more differential measurements (e.g. polarization) 
and/or with more extended kinematic coverage can allow us to disentangle 
different models and pin down the quarkonium production mechanism. 

With significantly improved statistics at high pT, RHIC experiments 
are now able to study the pT dependent J/ Raa [3 HOl IH] ■ Figure [s] 
shows the J/ip Raa results vs. centrality (left) and px (right). On the left 
plot, the latest STAR Au-|-Au results were divided into two px bins. The 
low pt (2-5 GeV/c) data show consistency with the PHENIX published 
data points. High px (> 5 GeV/c) data show systematically higher Raa 
compared to low px- In peripheral Au-|-Au collisions, the high px data are 
consistent with no suppression, while in central collisions, J/ip^s are still 
significantly suppressed, which may be due to the color-screening effect. On 
the right, the Raa values are plotted vs. px for two centrality bins and 
one can see the px dependent structure more clearly. Also plotted on these 
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Fig. 3. Left: Nuclear modification factor Raa of J/'0 vs. centrality. Measure- 
ments include PHENIX and STAR low results in Au+Au collisions, STAR high 
Pt results in Cu+Cu and Au+Au collisions. Curves shown in the plot are two 
model calculations which include both J/t/j dissociation and recombination pro- 
cesses. Right: Raa of J/V' vs. in Au+Au collisions from RHIC measurements 
compared to model calculations. 



two plots are theoretical model calculations which both include the J/V' 
dissociation in the QGP phase as well as the regeneration process [T2l ITS]. 
Both model calculations generally describe the data well. The system size 
and -pT dependence of J/V' Raa can be attributed to the formation time 
and/or leakage effects. 

There have been many interesting results from LHC experiments re- 
cently. ATLAS and CMS also observed large suppression of high -pT J I iii 
central Pb+Pb collisions [14^ [T5] , and the suppression seems to be stronger 
than that of RHIC high p^^ <//V''s, consistent with more suppression in a 
larger size system. ALICE reported low "pT J/'^ Raa at forward rapid- 
ity [1^. Comparing to the PHENIX result, the higher Raa values observed 
is consistent with more regeneration at LHC than at RHIC. Although cur- 
rent observations seem to qualitatively agree with expectations, but quan- 
titatively, we need further systematic studies on both RHIC and LHC to 
understand the quarkonium production mechanism as well as the medium 
properties. This also requires precision measurements in p[d)+A collisions 
to control the shadowing and the cold nuclear absorption effects. 

Another striking result reported by STAR is the J/tp elliptic flow (^2) 
measurement |17j . It has a significant improvement in term of precision 
compared the previous PHENIX measurement [18'. Figure |4] shows the 
J/V' V2 results from both PHENIX and STAR in 20-60% Au+Au colhsions 
and they are compared to various model predictions. Data show at pT > 2 
GeV/c, there is no sizable V2 for J/V', which disfavors coalescence production 
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Fig. 4. J/ip elliptic flow {V2) measurements from RHIC experiments in 20-60% 
Au+Au collisions compared with various model predictions. 



from thermalized charm quarks. 

One should note the valid centrality and pT regions for this statement. 
Based on model calculations [12] which generally reproduce the J/V' Raa 
and V2i the coalescence contribution in 20-60% Au+Au collision at RHIC 
is not dominant, and it mostly contributes in low (< 3 GeV/c) region. 
High pt region is still dominated by the initial J /ip production plus possible 
dissociation in the medium. To get insight of the clean charm quark V2 from 
the measurement of J/ ip at RHIC energy, one need to focus on the pT region 
below ~ 3 GeV/c and central collisions. This requires further improvement 
in the experimental precision. Another alternate way to learn the charm 
quark V2 will be to measure the charm hadron V2 as experimentally, and it 
will be achievable with great precision with help of the future silicon vertex 
detector. 

2.3. T measurements 

Study of bottomonium production at RHIC is a unique way to get insight 
of the originally proposed QGP color screen signature. The challenge for 
RHIC experiments is the statistics which makes this measurement likely 
to be a multi-year program. Both PHENIX and STAR have reported T 
(1S+2S+3S states if not specified) signals from p + p and Au+Au collisions 
since QM2009. The STAR pubhshed and PHENIX preliminary T cross 
sections in p + p collisions are consistent with pQCD model calculations 
and lie on the trend of energy dependence curve [IHl [20] . Recently, with 
significantly improved statistics in Au+Au collisions, STAR was able to 
measure the centrality dependence of T Raa |2I] • 

Figure [5] shows the T Raa vs. centrality. It shows a trend of T suppres- 
sion in central Au+Au collisions. Taking into account all current uncertain- 
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Fig. 5. T(1S+2S+3S) Raa vs. centrality from STAR compared to expectations 
from two naive T melting scenarios. 



ties from Au+Au and p + p reference, the measured value is ~ 3cj below 
unity. The red dotted line depicts a naive expectation assuming 2S and 3S 
states completely melt based on pQCD cross sections and PDG branching 
ratios. The blue dashed line depicts the expect value when all excited bot- 
tomonium states (including Xb) ah melt and only T(1S) survies [22]. The 
data now seem to favor the scenario that only T(1S) survies, but to learn 
quantitatively suppression levels of different bottomonium states requires 
significantly improved statistics in the future. 

With great mass resolution, CMS collaboration was managed to sep- 
arate (2S-I-3S) states from IS state and they reported the separate Raa 
values in minbias Pb+Pb collisions at 2.76 TeV [23]. When combining 
them together to compare with RHIC results, the expected Raa for central 
collisions should be less than 0.42, which is generally in line with the STAR 
result. However, one should note that the cold and hot medium effects may 
be significantly different between RHIC and LHC. 



3. Summary and Outlook 

In a short summary, RHIC heavy flavor program has achieved many 
significant results particularly in the past couple of years. It will remain as 
one of the focuses of RHIC heavy ion program in the upcoming RHIC II 
era. Both PHENIX and STAR are building significant detector subsystem 
upgrades now or in the coming years, e.g. PHENIX VTX and FVTX, STAR 
HFT and MTD upgrades etc. These are all aiming for precision measure- 
ments of both open heavy flavor and quarkonium production at RHIC II. 
With the LHC experiments ongoing, RHIC heavy flavor program will be 
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complementary and remain competitive in many aspects. Some of the mea- 
surements will be unique at RHIC, including: a) high precision open charm 
hadron measurements at low pT to address charm- medium interactions, b) 
bottomonium production measurements as bottomonia are expected to be 
clean at RHIC because of negligible contribution from regeneration, c) heavy 
quark correlation measurements as heavy quarks are expected to be back- 
to-back correlated in p + p collisions which allows clean interpretations for 
results in heavy ion collisions. These systematic study of heavy flavor mea- 
surements at RHIC and LHC will significantly improve our understanding 
of the sQGP matter by quantifying its physical properties with controlled 
accuracy. 

REFERENCES 

[1] Y. Zhang (for STAR Collaboration), QM 2011 proceedings, arXiv: 1106.6078. 

[2] M. Cacciari, P. Nason and R. Vogt, Phys. Rev. Lett. 95, 122001 (2005). 

[3] D. Acosta, et al. [CDF Collaboration], Phys. Rev. Lett. 91, 241804 (2003). 

[4] A. Dainese et al. [ALICE Collaboration], arXi v:1012.3 257, 

[5] B.I. Abelev et al. [STAR Collaboration], Phys. Rev. D 83, 052006 (2011). 

[6] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 97, 252002 (2006). 

[7] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 232301 (2007). 

[8] B.I. Abelev et al. [STAR Collaboration], Phys. Rev. C 80, 041902(R) (2009). 

[9] A. Adare et al. [PHENIX Collaboration], Phys. Rev. D 82, 012001 (2010). 
[10] Z. Tang (for STAR Collaboration), QM 2011 proceedings, arXiv: 1107.0532. 
[11] C. Powell (for STAR Collaboration), EPS 2011 proceedings, arXiv: 
1111.6944. 

[12] Y. Liu et al, Phys. Lett. B 678, 72 (2009). 
[13] X. Zhao and R. Rapp, Phys. Rev. C 82, 064905 (2010). 
[14] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 697, 294 (2011). 
[15] C. Silverstre (for CMS CoUaboration) , QM 2011 proceedings, arXiv: 
1108.5077. 

[16] P. Pillot (for ALICE Collaboration), QM 2011 proceedings, arXiv: 1108.3795. 
[17] H. Qiu (for STAR Collaboration), these proceedings. 
[18] A. Franz (for PHENIX Collaboration), J. Phys. G 35, 104002 (2008). 
[19] B.I. Abelev et al. [STAR Collaboration], Phys. Rev. D 82, 012004 (2010). 
[20] M. Leitch (for PHENIX Collaboration), Nucl. Phys. A 830, 27c (2009). 
[21] R. Reed (for STAR Collaboration), QM 2011 proceedings, arXiv: 1109.3891. 
[22] T. KoUegger, Ph.D. Thesis, Univ. of Frankfurt, 2005. 

[23] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 107, 052302 
(2011). 



